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Abstract. We address the estimation of the coupling constant of the 
£SJ Jaynes-Cummings Hamiltonian for a coupled qubit-oscillator system. 

! . We evaluate the quantum Fisher Information (QFI) for the system 

undergone the Jaynes-Cummings evolution, considering that the probe 
initial state is prepared in a Fock state for the oscillator and in a generic 
pure state for the qubit; we obtain that the QFI is exactly equal to the 
number of excitations present in the probe state. We then focus on the 
two subsystems, namely the qubit and the oscillator alone, deriving 
__ the two QFIs of the two reduced states, and comparing them with the 

previous result. Next we focus on feasible measurements on the sys- 
O | tem, and we find out that if population measurement on the qubit and 

1 1 1 Fock number measurement on the oscillator are performed together, 

the Cramer-Rao bound is saturated, that is the corresponding Fisher 
Information (FI) is always equal to the QFI. We compare also the per- 
formances of these measurements performed alone, that is when one of 
the two subsystem is ignored. We show that, when the qubit is pre- 
pared in either the ground or the excited state, the local measurements 
I are still optimal. 

> 

(Sj 1 Introduction 

OO 

^sO The Jaynes-Cummings (JC) model [JJ is one of the paradigmatic examples of hybrid 

systems, where a two- level system, modelled by a spin- 1/2, is coupled to a quantized 
y—{ mode of a harmonic oscillator. This is the typical situation in quantum electrodynam- 

t-H ics (QED) [2] where a single atom can be coupled to a cavity mode, for example in the 

t-H microwave [3|4] and in the optical [5] regime. The same model describes accurately 

other interesting physical systems, such as a single ion or a neutral atom in a trap [B], 
• i-H and the interaction of artificial atoms with resonators in circuit QED systems |7l8j . 

All these hybrid systems are considered essential for the development of quantum in- 
formation processing and in general of future quantum technologies |9|10j . The aim 
of quantum information is to characterize the peculiar properties of quantum systems 
and exploit them to perform tasks that would be not achievable in a classical context 
[IT]. For these purposes it is necessary to characterize the value of quantities that 
are not directly accessible either in principle or due to experimental impediments. 
This is the case of relevant quantities like phase [12113114] . entanglement |15ll6j or 
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temperature [17] , that cannot correspond to any quantum observable or also the cou- 
pling constant of different kinds of interactions [18 1 9120] . In this paper we apply the 
techniques of local quantum estimation theory (QET) |24l25l26l27j to the problem of 
estimating the coupling constant of the JC Hamiltonian. In this framework, the char- 
acterization of the estimation of the parameter is provided by the Fisher information 
(FI) which represents an infinitesimal distance between probability distributions and 
gives the ultimate precision attainable by an estimator via the Cramer-Rao theorem. 
Its quantum counterpart, the quantum Fisher information (QFI), is related to the de- 
gree of distinguishability of a quantum state from its neighbours [21122123] and gives 
the ultimate bound to the precision on the estimate allowed by quantum mechanics. 
In our estimation problem we consider the initial probe state prepared in a Fock state, 
as regards the oscillator and in a generic pure state, as regards the qubit. We address 
the overall estimation properties, by evaluating the QFI for the whole system under- 
gone the JC evolution. We also focus on the two subsystems alone, i.e. we consider 
the qubit subsystem obtained as the partial trace over the harmonic oscillator and 
viceversa, and we evaluate the QFI of the corresponding reduced states, in order to 
identify how much information about the parameter of interest is contained in each 
subsystem. Moreover we consider as possible feasible measurements on the coupled 
system, the measurement of the population of the excited state performed on the 
qubit and a measurement of the Fock number performed on the harmonic oscillator. 
We evaluate the FI for the collective measurement and observe that it allows one 
to achieve the ultimate bound on the precision given by the quantum Cramer-Rao 
bound. Finally we consider the FI for the same measurements performed on the two 
corresponding subsystems alone, that is when one of the two subsystems is ignored. 
We show that, when the qubit is prepared in either the ground or the excited state, 
both single local measurements are still optimal, that is the corresponding FIs are 
equal to the QFI of the whole systems. 

The paper is structured as follows: in Sec. [2] we describe the JC model and the uni- 
tary dynamics for the coupled system. In Sec. [3] we review some concepts of QET 
and illustrate the quantum Cramer-Rao bound focusing to the case of a pure-state 
unitary family. In Sec. [4] we show in details our results and finally, Sec. [5] gives some 
concluding remarks. 



2 The model 

The JC model describes the interaction between a single-mode bosonic field and a 
two-level system (qubit). Its dynamics is exactly solvable and the model has been 
widely investigated experimentally [5] . 

The bosonic field is described upon introducing an annihilation and a creation 
operator, satisfying [a,a^] = 1, and the corresponding Fock states {|^)}^L > *- e - tnc 
eigenstates of the number operator N = o)a, which provide a basis of the infinite- 
dimensional Hilbert space. In the rest of the manuscript we will adopt the quantum 
optics terminology calling photons the bosonic excitations, however the results pre- 
sented are also valid for the other physical settings described by the JC model. 
The two-level system (qubit) is characterized by a ground state \g) and an excited 
state |e) which are eigenstates of the Pauli operator a z and in turn form a basis for the 
two-dimensional Hilbert space of the qubit. Any pure state of the whole system com- 
posed by the bosonic field and the qubit can be written as = ^2j =g e S ra c j,n\j, n )' 
where we denote with \j, n) = \j) ® |n) the tensor product between the state of the 
qubit and the state of the field. 
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The JC Hamiltonian reads 

JiWq + ti-9 / + \ 

n = ^^a z + niLif a" a + -^-(cr + a + cj-a' ) (1) 

where the operators a± are the qubit ladder operators <j+ — |e)(<?| and er_ = |.g)(e|. 
The two operators a+a and o-a) correspond respectively to a transition from the 
lower level \g) to the upper level |e) together with the emission of a photon, and the 
transition from \g) to |e) together with the annihilation of a photon. It is clear that 
the interaction couples, for a given integer n, the states |e, n) and \g, n + 1). 

Upon choosing a suitable rotating frame and by considering the resonance con- 
dition u q = ujf, one rewrites the Hamiltonian % in the so called interaction picture 
as 

Hjc^Y^+ a + a - at) - (2) 
The corresponding evolution unitary operator for an interaction time r reads 

Ujc(n) = exp (- 1 -U jc t\ = exp(-WG) (3) 



where we defined 

G JC = and f2 = gr. (4) 



.at 



The parameter fl is the quantity of interest when we want to control the JC dynamics 
and in the following sections we will focus on its estimation properties. 

In our treatment we assume that at time t — the probe state is prepared in a 
pure state and no initial correlations between the qubit and the field are present, in 
formula 



e(o) = |!Pb)<*b| (5) 

with Itf'o) = ® l^/)- I n particular the qubit at time t = is prepared in a pure 
superposition of ground and excited states, 

9 6 

\ip q ) = cos-|e) + sin (6) 

while the bosonic field is prepared in a Fock state 1-0/) = \n). Notice that preparation 
of Fock states have been proposed theoretically and proved experimentally both in 
cavity QED, ion trapping and circuit-QED systems [2I29I30I31I32I33I34I35] . Given 
this probe preparation, the evolution of the system reads 

g(O) = Uj C (n)g(0)u] c (f2). (7) 

and upon tracing the evolved state over the bosonic field or the qubit degrees of 
freedom, we obtain the states 

eq (n) = Tr f [Ujc{f2)Q(p)rt JC (!2)] (8) 
Q f (S2) = Tr q [U JC (O)Q(0)Ul c (n)} (9) 

which describe respectively the qubit state and the harmonic oscillator subsystem 
state at time t. 
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In details, the reduced qubit density operator reads in the basis {|e), \g)} as 

*<°> -(££)■ <10) 

where 

2 2 I ^ I A 2 ^ 2 f ^ 

g gg =cos -sin I— Vn + ll+sin -cos I t^V™ 

Qee — 1 Qgg 

Beg = ^ sin 6 1 cos (J^VnJ cos (J^^ n + = Qge- (H) 

On the other hand, the reduced density operator for the bosonic field is diagonal in 
the Fock basis: 



g f (f2) =p n _ 1 \n- l)(n - 1| +p n \n)(n\ + p n+1 \n + 1) (n + 1|, (12) 



where 



. 2 6> /l-cos(l2v/n)\ 
Pn-i = sm 2 - j 



1 / 6* 6* 

p„ = - ( 1 + cos2 2 cos + l) + sin 2 - cos 



2 / 1 - COS l) 



Pn+1 = COS 2 ( 2 J ' ^ 



3 Local quantum estimation theory 

In this section we review the basic concepts of local quantum estimation theory that 
will be applied later to the qubit-oscillator system. 

An estimation problem consists into choosing a measurement somehow related 
with the parameter of interest and then defining an estimator, i.e. a function from 
the set of the measurement outcomes to the parameter space, in order to infer the 
value of the quantity that we want to estimate. Classically, given the conditional 
probability p(j\0) of measuring the outcome j when the value to be estimated is £2, 
optimal estimators are those saturating the Cramer-Rao inequality, which establishes 
that the variance Var(J?) of any unbiased estimator is lower bounded by 

Var ^) * MiW (14) 

where M is the number of measurements of the sample and F{ff) the Fisher infor- 
mation (FI) 

F{Q) = Y,P{m[do\np{j\Q)] 2 . (15) 

3 

In quantum mechanics, according to the Born rule one has p(j\f2) = Tr[gnIIj], where 
go is a family of quantum states which depend on the parameter SI and the operators 
{II j} are the elements of the probability operator- valued measure (POVM) describing 



Will be inserted by the editor 



5 



the quantum measurement. By defining the symmetric logarithmic derivative (SLD) 
operator Lq by means of the following equation 

dgn LqQq + qqLq 

M2= 2 ' (16) 



the classical Fisher Information in Eq. (pL5|) can be rewritten as 

Tr[f> fi i7 



F ^> ~ — tQi — (17) 



which establishes the maximum precision for the estimation of the parameter fl for 
a fixed measurement {II j}. Moreover, maximizing the FI over all the possible mea- 
surements, one can show that 

F(f2) < H{fi) = Tc[ 8n L%]. (18) 

The quantity H(Q) is called quantum Fisher information (QFI) and define the cor- 
responding quantum Cramer-Rao bound 

Var(l2) > 1 —, (19) 

which gives the ultimate limit to the precision allowed from quantum mechanics for 
the estimation of the parameter fl labelling a given quantum statistical model qq. 

When the set of quantum states qq is given in a diagonal form qq = ^2 k Xk\4>k) (<fik\, 
the QFI can be evaluated as 

H(f2) = £ + 2^ {X ^f\(Md^ k )\ 2 . (20) 

A particular case is given when the family of qq is a unitary pure-state family, 

Qn = U n \W )(%\uj } U n = expHG/2) (21) 

where G is the generator of the transformation and the initial pure probe state. 
In this case one can prove that the QFI is independent on the actual value of the 
parameter 12 and it is proportional to the fluctuations of the generator on the probe 
state, i.e. 

H(f2) = 4 (A 2 G) 

= 4((^o|G 2 |^o)-(^o|G|^o) 2 ) (22) 

One can thus rewrite the quantum Cramer- Rao bound and read it as a Heisenberg-like 
uncertainty relation 

Var(tf) (A 2 G) > ~ (23) 



4 Results 



In this section we report the results for the estimation of the parameter Q in the JC 
model. We first derive the ultimate limit to the estimation by evaluating the QFI 
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for the whole coupled system composed by the qubit and the harmonic oscillator. 
We then consider the two subsystems alone, by performing a partial trace either on 
the bosonic field or the qubit, and evaluate the QFI of the reduced states Q g (f2) and 
Qf(£2) given respectively in Eqs. (10) and (12) . We also consider the performances 
of two feasible measurements: the photon number measurement on the bosonic field 
and the population measurement on the qubit. We first evaluate the FI for the col- 
lective measurement on the total system and finally, we consider the FI for the same 
measurements performed on the two subsystems alone, namely the FI for population 
measurements on the qubit subsystem and the FI for photon number measurements 
on the oscillator subsystem alone. 



4.1 Quantum Fisher Information 

In the following we evaluate the QFI for the pure state g(f2) obtained upon choosing as 
a probe state |#q) = \ip q ) <S> \n), where the state of the qubit is a general superposition 
\ip q ) as in Eq. (T6J) and the field is prepared into a Fock state |n). 
We consider the evolution of the whole system and assume that both the degrees of 
freedom of the qubit and of the field are accessible and measurable, i.e. our statistical 
model is a pure state unitary fam ily and the corresponding QFI for the parameter 



fi can be evaluated as in Eq. (22) with the generator Gjc = + aV_)/2. The 
calculation leads to 

Q 

H(Q) = n + cos 2 -. (24) 

Note that the QFI H(fl) has the minimum value for 8 = ir, where H(f2) = n, and 
the maximum is H{S1) = n+ 1 for 6 = 0. Thus the optimal preparation for the qubit 
state which gives the maximum value of the QFI, corresponds to the excited state 
\ip q ) = \e). In particular one can observe that the QFI H(Q) of the total system is 
equal to the total number of excitations of the probe state, 

H{Q) = (&o\E\<P } where E = a) a + cr+cr_, (25) 

showing that the more excitations we have, the more precise will be the estimation 
of the JC coupling constant. 

We now consider the case where the degrees of freedom of one of the two subsystems 
are not accessible. The reduced states of the qubit and the bosonic field are given 



respectively in Eqs. (10) and (12) and the corresponding QFIs are denoted as H q {Q) 



and Hf(Q). We evaluate theniBy means of Eq. (20) and plot their behaviour in Fig. 
[T]as a function of 9 and for different values of the parameter SI. We observe that the 
bosonic field alone contains more information about the parameter SI than the qubit 
subsystem, that is 

H f {f2)> H q {Q). (26) 

Moreover, we also observe that if the qubit is prepared in either the ground or excited 
state, the values of the QFI of the total system H(S1) coincides to the QFIs of the 
two reduced subsystems. In formula we have that, for all the possible values of J?, 

H q (O)[e = 0] =H f (Q)[9 = 0] =ra+l (27) 
H q (f2)[6 = tt] = H f (Q)[6 = 7r] = n. (28) 

This result shows that, for particular choices of the input state, a measure performed 
on one of the subsystems, either the qubit or the oscillator, can attain the ultimate 
precision on the estimate of the parameter SI given by the QFI H(fl). 
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H(H) 




Fig. 1. QFIs of the two subsystems as a function of 6 with n — 3 and for different values 
of the parameter Q (blue lines, Q = f.O; red lines, Q = 1.5). The dashed lines denote the 
QFI Hf(H) of the bosonic field, while the solid lines denote the QFI H q (Q) of the qubit 
subsystem. The black solid line on top, represents the QFI H(O) of the qubit-oscillator total 
system. 



4.2 Fisher information for Fock and population measurements 

In order to address the performances of feasible measurements for the estimation of 
the parameter J?, we consider the FI for the photon number measurement on the 
bosonic field and the population measurement on the qubit system. If we consider 
a collective measurement performed on both subsystems, the non-zero conditional 
probabilities that we have to take into account are 

(29) 
(30) 
(31) 
(32) 

where p(j, n\(7) denotes the conditional probability of obtaining the qubit in the state 
j and the bosonic field with n excitations, when the parameter has the value Q. The 
corresponding Fisher information is evaluated by means of Eq. ( 15 ) giving 

Q 

F(f2) = H{fi) = n + cos 2 -. (33) 

Therefore the Fisher information F{f2) for the Fock and population measurements 
saturates the quantum Cramer-Rao bound i.e. the measurement considered is always 
optimal. 

We now study the two measurements separately, i.e. we evaluate the Fisher infor- 
mation for population measurements on the qubit F q ({2) and the Fisher information 
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F(0)/H(n) 




Fig. 2. (From top to bottom) Black solid line: QFI H(Q) and FI F(Q) for collective 
Fock and population measurement. Dashed blue line: QFI Hf(Q) and FI Fj{Q) for Fock 
measurement on the harmonic oscillator subsystem alone. Dashed red line: QFI H q (Q). 
Dotted red line: FI F q (f}) for population measurement on the qubit subsystem alone. The 
parameters considered are fl = 1.25 and n = 3. 



for the photon number measurements on the harmonic oscillator Ff(f2). In particular 
we want to compare them with the corresponding QFIs for the two subsystems H q {Q) 
and Hj(f2). The population measurement on the qubit corresponds to measure the 
Pauli operator a z , whose eigenstates are indeed {|e), \g)}. The conditional probabil- 
ities for the qubit to be found in the ground or the excited state are given by the 
diagonal elements of the reduced state g q {fi), that is 

2 9 n I il i \ 9 ^ I — \ 

V\9\U) = Q gg = cos -sin I — Vn + 1 I + sin -cos I — V™ I (34) 

p(e\f2) = Q ee = 1 - Q gg . (35) 

On the other hand, the non-zero conditional probabilites for the number measurement 
on the bosonic fields are 

. , ,9/1- cos (flJn) 

p(n - l\fi) = p n _j = sin 2 - f 

1 / 6 9 

p(n\fi) — p n = — ( 1 + cos 2 - cos (i?V« + l) + sin 2 - cos (^?^ 

9 9 ( 1 - cos (f2y/n + l) \ 
p(n + 1| S2) = p n+1 = cos 2 - 1 (36) 



The corresponding FIs, F q (Q) and Ff(f2), are evaluated by means of Eq. (15) and 
are plotted in Fig. [2j Since the field subsystem state Qf{fi) is diagonal in the Fock 
number basis, we have that the Fock number measurement is always optimal, i.e. the 
Fisher information is equal to the quantum Fisher information of the state Qf{fi) 

F f (n) = H f ({2) v n,e,n. (37) 



On the other hand, the population measurement is in general not optimal but it 
saturates the subsystem quantum Cramer- Rao bound when 9 = or 9 = 7r, that 
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is when the qubit is prepared respectively in the excited state \e) or in the ground 
state \g). In these special cases all the quantities we have evaluated so far are equal, 
that is the FIs for the measurements on the two subsystems are equal not only to the 
corresponding QFIs of the reduced states Hf(fl) and H q (S7), but also to the QFI of 
the total system H(Q). In formula we have that 



This interesting feature suggests that for these particular choices of the qubit prepara- 
tion, in order to attain the ultimate limit posed by the quantum Cramer-Rao bound, 
we can decide to measure only one of the two subsystems neglecting the other degrees 
of freedom. 



5 Conclusions 



In this paper we addressed the estimation of the coupling constant of the Jaynes- 
Cummings Hamiltonian considering as an input probe state the field prepared into 
a Fock state |n) and the qubit prepared in a generic superposition of excited and 
ground state. 

We derived the QFI for the whole system, i.e. the qubit-oscillator system evolved un- 
der the JC Hamiltonian and observed that it is equal to the number of excitations of 
the input probe state and independent on the value of the parameter to be estimated. 
We then derived the QFI for the two subsystems, i.e. the QFI of the qubit and the 
oscillator state. In this case we have found that in principle the parameter is better 
estimated when the qubit subsystem is ignored, rather than the bosonic field. 
We then considered a feasible detection scheme where a measurement of the popula- 
tion of the excited state is performed on the qubit and a measurement of the Fock 
number is performed on the harmonic oscillator. The FI for the corresponding collec- 
tive measurement turned out to be equal to the corresponding QFI, saturating the 
quantum Cramer-Rao bound and providing the optimal estimate of the parameter. 
Finally we considered the FI for the same measurements performed separately on the 
two subsystems alone, namely by ignoring the degrees of freedom of one subsystem. 
The surprising result is that, if the qubit is prepared either in the ground or in the ex- 
cited state, both the measurements performed on the single subsystems, and ignoring 
the remaining one, are optimal. This result is relevant in many practical situations 
where one of the two subsystems is not experimentally accessible. Moreover since the 
bound obtained does not depend on the value of the parameter, it is not necessary 
to tune the measurement by means of two-step or adaptive estimation strategies in 
order to attain the optimality. 

As we stressed above, our study shows the existence of a link between the estimation 
of the JC coupling constant and the amount of excitations present in the probe state. 
Since the preparation of a high-number Fock state is still experimentally challenging, 
as a future outlook, it would be interesting to study different preparations of the 
probe states, for example by considering the bosonic field in a coherent state with 
a high number of photons. It will be relevant in this case, also to understand if the 
measurements that in this case has been proved to be optimal, remain optimal with 
such a different preparation of the probe. 




(38) 
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